Variational Monte Carlo approach to the two-dimensional Kondo lattice model 
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We study the phase diagram of the Kondo-lattice model with nearest-neighbor hopping in the 
square lattice by means of the variational Monte Carlo technique. Specifically, we analyze a wide 
class of variational wave functions that allow magnetic and superconducting order parameters, so 
to assess the possibility that superconductivity might emerge close to the magnetic instability, as 
often observed in heavy fermion systems. Indeed, we do find evidence of d-wave superconductivity 
in the paramagnetic sector, i.e., when magnetic order is not allowed in the variational wave function. 
However, when magnetism is allowed, it completely covers the superconducting region, which thus 
disappears from the phase diagram. 

PACS numbers: 71.10.Hf, 71.27. +a, 75.20.Hr, 75.30.Mb 



I. INTRODUCTION 

The Kondo lattice model (KLM) describes localized 
magnetic moments that interact with a single band of 
itinerant electrons via an antiferromagnetic exchange 
coupling J. This model has been introduced long ago 
by Doniachji and, since then, has been widely invoked to 
describe the physics of heavy-fermion materials?- 2 Three 
ingredients define the KLM: the bandwidth W of the 
conduction electrons (usually denoted by c electrons), 
their density n c , and the Kondo exchange J that con- 
trols the coupling between the local moments and the 
spin density of c-electrons. The localized moments (de- 
noted by / electrons) are anchored to the sites of a reg- 
ular lattice. The phase diagram of the KLM depends 
in a non-trivial way on J/W and the electron density 
n c . In one spatial dimension, the KLM has been inten- 
sively studied and shows three distinct phases 4 In the 
compensated case, one conduction electron per localized 
spin, i.e., n c = 1, it is a so-called Kondo insulator with 
a charge as well as a spin gap without any magnetic or- 
der, a kind of spin-liquid insulator. For n c ^ 1, it is 
either a paramagnetic metal for low J/W, or a ferromag- 
netic metal for larger J/W. In higher dimensions, where 
the SU(2) spin symmetry can be spontaneously broken, 
the phase diagram is expected to enrich and display a 
critical point separating a paramagnetic heavy-fermion 
metal from a different metallic phase with long-range 
magnetic order Indeed, in analogy with the Kondo 
effect that occurs in the case of a single magnetic impu- 
rity, the conduction electrons may screen the localized 
moments, thus forming a global singlet state. However, 
such a Kondo screening is thwarted by the tendency to 
magnetic ordering of the localized moments. In fact, the 
latter ones interact mutually via the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) mechanism, through an effec- 
tive exchange J e ff(q) oc — J 2 3?ex(<z,w = 0) mediated 
by the magnetic polarization xili w ) °f the conduction 
band. The competition between Kondo screening (that 
favors a paramagnetic ground state) and RKKY interac- 
tion (that favors magnetically ordered states) is the heart 



of a frustrating behavior, which may lead to genuine 
quantum phase transitions. 1 Furthermore, other physical 
processes may profit from the balanced competition be- 
tween Kondo screening and magnetic ordering nearby the 
transition, and make novel phases to intrude, most no- 
tably superconductivity. Indeed, the discovery of super- 
conductivity in CeCu2Si2^ and subsequently in several 
other heavy-fermion materials, unveiled the rich variety 
of phenomena of these strongly-correlated systems. It 
is widely believed that, in heavy-fermion compounds, su- 
perconductivity is not the conventional phonon-mediated 
one, but most likely it is caused by antiferromagnetic spin 
fluctuations^ 

Recent experiments of the Hall coefficient in the heavy- 
fermion material YbRh2Si2j£ have provided new intrigu- 
ing elements that renewed interest in the phase diagram 
of the KLM. In particular, the rapid jump of the Hall 
coefficient^— is suggestive of a sudden change of the 
Fermi surface topology at, or close to, the quantum phase 
transition between the non-magnetic metal and the mag- 
netic one. In the conventional view of heavy-fermions, 
the localized spins are promoted in the conduction band 
through the Kondo effect. A mass-enhanced Fermi liq- 
uid is thus settled down, with a "large" Fermi surface 
that includes the conduction as well as the localized elec- 
trons. In the simplest scenario of a spin-density-wave 
quantum phase transition ; 14 i 15 magnetic ordering is in- 
deed expected to modify the topology of the Fermi sur- 
face by appearance of the magnetic Bragg reflections. 
In addition, the Fermi surface reconstruction across the 
magnetic field induced transition in YbRfi2Si2 can also 
be simply interpreted as a Zeeman driven Lifshitz tran- 
sition of the heavy quasiparticles. 16 However, an alterna- 
tive scenario is possible in which the Kondo effect dies 
out at the transition point ^ hence the local moments 
suddenly stop contributing to the volume of the Fermi 
surface, which then counts only the number of c elec- 
trons. In the language of the Anderson lattice model, 
which maps for strong repulsion onto the KLM, the death 
of the Kondo effect would translate into the Mott local- 
ization of the / electrons, whose magnetic ordering would 
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then be only a by-product rather then the driving source. 
Indeed, there are suggestions that the Kondo-breakdown 
and the onset of magnetic order in the Anderson lat- 
tice model are distinct phenomena, which may occur at 
different points of the phase diagrami 18 i 19 This scenario 
has been indirectly supported by variational Monte Carlo 
calculations 20 and by the Gutzwiller approximation 2 ^ in 
the KLM on a square lattice. Indeed, these works found 
evidence of two distinct phase transitions, one given by 
the continuous appearance of long-range magnetic order 
and another one related to an abrupt topological change 
of the Fermi surface. The same outcome of two distinct 
transitions has been later observed also within the dy- 
namical cluster approximation (DC A). 22 In spite of all 
efforts, such an interesting issue like remains open. 

In this paper, we investigate the KLM in two dimen- 
sions paying attention not only to magnetism, but also to 
the possible emergence of superconductivity in its prox- 
imity. We use both mean-field and variational Monte 
Carlo approaches. As far as the former one is concerned, 
an analytical treatment of the long-range antiferromag- 
netic phase is possible only at compensation, n c = 1, 
where calculations have been already performed, 23 Here, 
we go beyond the results of Ref. and consider also the 
uncompensated regime, n c < 1, by solving numerically 
the Hartree-Fock equations. Furthermore, we general- 
ize the previous calculations based upon the variational 
Monte Carlo approach— or the Gutzwiller approxima- 
tion 2 ^ and introduce additional correlations inside the 
trial wave function, among which superconducting ones. 
In particular, our variational calculations show that su- 
perconductivity is indeed stabilized in the paramagnetic 
sector in a region of the phase diagram close to n c = 1 
and not too large J/W. However, when magnetism is 
allowed in the variational wave function, an antifcrro- 
magnetic phase emerges and completely covers the super- 
conducting dome. Therefore, at least within the varia- 
tional approach and in our model where the only source of 
magnetic frustration is deviation from the compensated 
regime, we do not find any superconducting region at 
the border between antiferromagnetic and heavy-fermion 
metal. 

The paper is organized as follow: in section [Hi we in- 
troduce the microscopic model and discuss the methods 
that we use; in section UTTl we present our numerical re- 
sults for both mean-field and variational approximations; 
and, finally, in section ITVl we draw our conclusions. 



II. MODEL AND METHODS 

The KLM model on the two-dimensional square lattice 
is defined by: 

n = -tJ2 c Ls> + h - c - + J E s * ■ s * w 

where denotes nearest-neighbor sites i and j and 

c l a ( c i a) creates (destroys) an itinerant electron at site 



i with spin er; s, = (sf,sf,sf) is the spin operator for 
the c-electrons, i.e., sf = 1/2 £ CT Cr , c\ a T* a ,c i <yl , r a be- 
ing the Pauli matrices. Similarly, S, = (Sf,S?,S£) is 
the spin operator for the localized / electrons, Sf = 
1/2 V , f] „t"„, f„- „/. By constraint there is one f elec- 
tron per each site. The exchange coupling is antiferro- 
magnetic, i.e., J > 0, and we shall take all energies mea- 
sured in units of t. 



A. Mean-field approach 

The simplest approach to the KLM of Eq. ([1} is the 
mean-field approximation, in which the spin-spin interac- 
tion is decoupled to bring about a non-interacting Hamil- 
tonian. We shall implement Hartree-Fock by assuming 
finite the following average values, to be determined self- 
consistently: 



m f = (-!)«' (Sf). 



(2) 

(3) 
(4) 



V denotes the (site-independent) c—f hybridization, re- 
sponsible in mean-field for the creation of the Kondo sin- 
glet; m c and m/ are the staggered magnetizations of con- 
duction and localized electrons, respectively, which have 
opposite sign because of the antiferromagnetic exchange 
J. In addition, a Lagrange multiplier /i/ must be in- 
cluded to enforce (on average) the /-orbital occupancy, 
i.e., rif = 1. 

In momentum space, the antiferromagnetic mean-field 
Hamiltonian can be cast in a 4 x 4 matrix form: 
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(5) 



where the sum over k is restricted to the reduced (mag- 
netic) Brillouin zone. In order to compute the total en- 
ergy, the constant term (3JV 2 /2 + Jm c m,f + [if)N must 
be added (N being the number of sites). 

The paramagnetic state is found by imposing m c = 
m f — and corresponds to the 2x2 matrix Hamiltonian: 



H 



PM 



E 

k.a 



[ c k.a Jk,i 



V - M/ 



fk. 



(G) 



The self-consistency conditions Eqs. ©, Q, and (0| 
are solved numerically on finite size systems with N 
sites, number that we scale to get reliable estimates in 
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FIG. 1. (Color on-line) Quasi-particle band structure for var- 
ious phases: the AFs (a) has flat / bands, due to the full de- 
coupling between / and c electrons and a small Fermi surface; 
the antiferromagnetic phases with electron-like and hole-like 
Fermi surfaces are denoted by AFe (b) and AFh (c), respec- 
tively; finally, PM (d) is a paramagnetic (metallic) phase. The 
dashed red line indicates the chemical potential. 



the thermodynamic limit. We mention that an ana- 
lytic solution of the problem is possible only in the com- 
pensated regime, while in general numerical calculations 
are needed. In practice, we numerically diagonalize the 
4x4 matrix for all k points independently and then fill 
the bands with the lowest orbitals; the mean-field pa- 
rameters are numerically calculated and the procedure 
is iterated until convergence is reached. At the mean- 
field level, a superconducting singlet order parameter 
So- a ( C I a f I - a ) i s n °t independent from the hybridiza- 
tion, because of the charge-isospin SU(2) symmetry dis- 
played by the / electrons. 

Within the AF state, we find three different cases, de- 
pending on the magnitude of variational parameters. In 
the following, we will adopt the notations of Ref. I20I . 
Whenever the hybridization parameter vanishes, i.e., 
V = 0, the localized electrons decouple from the con- 
ducting ones and they do not contribute to the volume 
enclosed by the Fermi surface, in this case we have an 
antiferromagnetic state with a "small" Fermi surface (de- 
noted by AFs). By adding a small hybridization to the 
AFs, we end up with a state which has an electron-like 
Fermi surface, the so-called AFe. Finally, in the case 
where the hybridization V is large and the magnetic or- 
der parameter small, we have a hole-like Fermi surface, 
the so-called AFh. Here, the / electrons participate to 
the total volume enclosed by the Fermi surface, which is 
therefore "large" . A qualitative picture of all these cases, 
together with the simple paramagnetic state is depicted 
in Fig. HJ 



B. Variational wave functions 

In order to go beyond the mean-field approximation, 
we consider correlated variational wave functions, in 
which the constraint of one / electron per site is im- 
posed exactly via a Gutzwiller projector. This is achieved 
through the projected variational wave function: 

I*) =7>/|*mf), (7) 

where Vf is the projector which enforces single occupa- 
tion of / orbitals on each site. Here, |\&mf) is an un- 
corrected wave function defined as the ground state of a 
non-interacting variational Hamiltonian Kmf that may 
contain, in addition to the mean-field parameters of the 
previous section iHAl also direct /— / hopping as well as 
superconducting terms: 

xfj = (flU,}, (8) 

A SH c !t c L + c iT c a>' (io) 

in s-wave or o?-wave configurations. An on-site c—c pair- 
ing has been also considered. 

Therefore, in the following, we shall consider four kind 
of uncorrelated variational wave functions: (1) paramag- 
netic, (2) antiferromagnetic, (3) superconducting, and, 
finally, (4) with coexisting antiferromagnetism and su- 
perconductivity. The variational parameters of the non- 
interacting Hamiltonian are determined so as to mini- 
mize the total energy. Because of the presence of the 
Gutzwiller projector "P/ we have to use a variational 
Monte Carlo technique 24 to compute the total energy. 
In practice, we minimize the variational energy for all 
the previous states as a function of the exchange cou- 
pling J and the electron density n c . Calculations have 
been performed on clusters with 64, 100, 144, and 256 
sites. Suitable boundary conditions have been chosen to 
obtain close-shell configurations in I^mf)- 

III. RESULTS 

Here we present our numerical results on the KLM, 
first within the mean-field approximation and then by 
the variational Monte Carlo approach. 

A. Mean-field results 

The mean-field phase diagram, as a function of J and 
the electron density n c , is reported in Fig. [5] (for a direct 
comparison, we report the variational Monte Carlo phase 
diagram in Fig. [3]). For n c > 0.81 we find two distinct 
phase transitions. When J is small, the ground state has 
antiferromagnetic long-range order and displays a small 
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FIG. 2. (Color on-line) Mean- field phase diagram. AFs 
and AFh indicate antiferromagnetic phases with small and 
hole-like Fermi surfaces, respectively. PM indicates a para- 
magnetic (metallic) phase. Continuous (dashed) lines mark 
second-order (first-order) phase transitions. 




FIG. 3. (Color on-line) Variational Monte Carlo phase dia- 
gram. The labels are the same as in Fig. [2] We note that, 
although qualitatively similar, the variational and mean-field 
phase diagrams are quantitatively quite different. 
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FIG. 4. (Color on-line) Mean-field order parameters (hy- 
bridization V and total magnetization m z = m/ — m c ) as a 
function of the Kondo exchange for two different densities of 
the conducting electrons, n c — 0.62 (left panel) and n c = 0.93 
(right panel). 



first order. In Fig. [U we show the behavior of the hy- 
bridization V and the total magnetization m z — m,f—m c . 
Eventually, by further increasing the local exchange, the 
Kondo mechanism prevails and the system becomes a 
paramagnetic metal where conduction electrons screen 
the local moments. The transition between the AFh 
phase and the paramagnetic metal is second order, with 
the magnetization that goes continuously to zero, see 
Fig. [4j Moreover, the hybridization is continuous though 
the transition and the topology of the Fermi surface does 
not change. 

For smaller values of the conduction electron density, 
i.e., n c < 0.81, the AFh state cannot be stabilized and 
there are only two phases: the AFs for small Kondo ex- 
change and the paramagnetic metal for large ones. The 
phase transition between them is first order: both the 
antiferromagnetic order parameter and the hybridization 
change abruptly from zero to a finite value, see Fig. 2J In 
this case, also the topology of the Fermi surface changes 
across the transition. 



Fermi surface, namely we obtain the AFs state. Here, the 
local / electrons are totally decoupled from the conduct- 
ing ones (the mean-field equations are solved by V = 0) 
and do not contribute to the Fermi surface. This regime 
is dominated by the RKKY interaction that generates a 
magnetic pattern in the localized spins, and consequently 
also in the conducting ones: the magnetization of / elec- 
trons is saturated, i.e., mj = 0.5, while m c is a smooth 
function, slightly increasing with J. 

By increasing J, the Kondo mechanism becomes com- 
petitive with the RKKY interaction and we enter into 
another antiferromagnetic phase, where c and / electrons 
are hybridized. Here, there is a hole-like Fermi surface 
and, therefore, the phase is AFh. The hybridization V 
has a finite jump at the transition, which is, therefore, 



B. Variational Monte Carlo results 

Now, we turn to the variational Monte Carlo results, 
summarized in the phase diagram of Fig. [3] (to be com- 
pared with the mean- field one, see Fig. [2]). Within the 
variational Monte Carlo, which improves substantially 
the Hartree-Fock energy, the c— / hybridization parame- 
ter V of the no n- interacting auxiliary Hamiltonian T-Lmf, 
is finite throughout the phase diagram^ It follows that 
the zero-temperature variational Fermi surface always in- 
cludes both c and / electrons. However, the optimized 
V is tiny in the AFe phase and steps up discontinuously 
entering the AFh or PM phases. Therefore, a very small 
temperature can wash away the effects of V in the AFe 
phase (but not in the AFh and PM ones) thus better 
highlighting the differences between the phases. For this 
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reason, we decided to follow Ref. |2l| and calculate, in the 
Brillouin zone, the emission spectrum A(k) of %mf at 
the chemical potential broadened with a low but finite 
temperature T: 



A(k) 



deA{k, 



dm 

' 8e ' 



where A(k, e) is: 



A(k,e) = K™kJ0)| 2 <5(e - E n + E ), 



(12) 



(13) 



n>0 



where \n) are the unprojected (mean-field) states, with 
energies E n . 

We start our analysis by considering the paramagnetic 
sector, which is reacher than the one obtained within 
mean-field approximation, and can shed some light by 
disentangling Kondo effect from long-range magnetism. 
The paramagnetic phase diagram of the KLM, allowing 
for superconductivity, is shown in Fig. [5] We find that, 
although (on-site or extended) s-wave pairing is never 
stabilized, a sizable ci-wave pairing is obtained in a wide 
range of parameters, namely for J < 1.5 and n c > 0.65, 
and brings a non-negligible energy gain with respect to 
a normal phase. The condensation energy is reported in 
Fig. [B] for three values of n c . For J < 0.1, the pairing 
correlations of the unprojected state become very small, 
implying a tiny energy gain with respect to the nor- 
mal state. We emphasize that superconductivity emerges 
only thanks to the electronic correlations brought by the 
Gutzwiller projector Vf, since pairing does not arise at 
the mean-field level. A finite ci-wave pairing is thus gener- 
ated by the antiferromagnetic c— / exchange, suggestive 
of similarities to analogous results found in t— J models 
for cuprate superconductors ! 25 ' 26 Indeed, as evident from 
Fig- El the condensation energy has a bell-like shape, with 
maximum at some intermediate values of J and n c . 

We mention that a very recent single-site dynamical 
mean-field theory (DMFT) calculation in the paramag- 
netic sector— finds evidence of s-wave superconductivity, 
whose maximum strength is reached, for a semicircular 
density of states, when J ~ 1.6 (translated in our units 
in which the bandwidth is 8t) and n c ~ 0.86, which we 
could not reproduce by our simple variational wave func- 
tion. 

In Fig. [7J we plot A(k) for four different values of J, 
two well inside the superconducting region and two across 
the transition to the normal phase. Since the transition is 
continuous, the T > Fermi surface continuously change 
from electron-like to hole-like, see Fig. [7] A large spectral 
weight along the zone diagonals in the superconducting 
phase is observed whenever sizable pairing correlations 
are present, because of <i-wave symmetry. 

When we leave the paramagnetic sector and allow for 
antiferromagnetism, the latter prevails over superconduc- 
tivity, which therefore disappear from the actual phase 
diagram, see Fig. [3] In other words, the energy gain 
of antiferromagnetism always overcomes that of super- 
conductivity, see Fig. [SI ruling out the possibility of a 
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FIG. 5. (Color on-line) Variational phase diagram in the para- 
magnetic sector. PM and SC denote the paramagnetic metal 
and the d-wave superconducting state, respectively. The tran- 
sition between these two phases is always continuous. 
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FIG. 6. (Color on-line) Energy (per site) difference between 
the superconducting state and the metallic one as a function 
of J for different values of the c-electron density. 



ground state with superconductivity and no magnetic or- 
der. This occurs at least in the bipartite nearest-neighbor 
hopping model that we have considered, where the only 
source of frustration is the conduction electron density 
n c lower than half-filling. We also investigated possible 
coexistence between antiferromagnetism and ci-wave su- 
perconductivity, which we indeed found but only in the 
AFe region. However, we believe this result is only a 
finite size effect, since the energy gain by allowing d- 
wave pairing on top of magnetism is tiny (at maximum, 
AE ~ 10~ 4 i) and, in addition, the size scaling of the 
actual order parameter (after Gutzwiller projection) sug- 
gests a vanishing value in the thermodynamic limit. We 
observe that the region of stability of the AFe phase is 
reduced substantially with respect to the corresponding 
AFs found at the mean-field level, compare Fig. [3] with 
Fig. [21 showing that the variational wave function can 
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FIG. 7. (Color on-line) Emission spectrum A(k), see Eq. (fT2jl . 
broadened with a temperature T = 0.01 for n c — 0.91. The 
top panels correspond to J = 0.2 (left) and J = 0.8 (right), 
the bottom one to J = 1.2 (left) and J = 1.3 (right). The 
lower three values of J are inside the d-wave superconducting 
dome, while J — 1.3 is already in the metallic phase. Note 
the change of topology as J increases. 



FIG. 9. (Color on-line) Emission spectrum A(k), see Eq. HU, 
broadened with a temperature T = 0.01. Top left panel: n c = 
0.9375 and J = 0.3, inside the AFe phase. Top right panel: 
n c — 0.9375 and J = 1.3, inside the AFh phase. Bottom left 
panel: n c = 0.75 and J = 0.2, inside the AFe phase. Bottom 
right panel: n c = 0.75 and J = 1.2, inside the PM phase. 
Note the shadow bands for the antiferromagnetic cases. 




FIG. 8. (Color on-line) Energy (per site) difference between 
the antiferromagnetic, superconducting, and normal states as 
a function of J for n c = 0.91. 



deal with Kondo screening better than mean held. 

In Fig. [9] we draw A(k) for different values of n c and 
J. The left panels are inside the AFe phase, and show a 
spectral distribution at the chemical potential that cor- 
responds to a small, electron-like, Fermi surface. On the 
contrary, the right panels (the top one inside the AFh 
phase and the bottom one in the PM region) indicate a 
larger Fermi surface that still contains / electrons at that 
value of temperature T. We note the signals of shadow 
bands in the antiferromagnetic A(k) of the top panels 
and left bottom one. 



We finally mention that the phase diagram of Fig. [3] 
agrees pretty well with that obtained by Watanabe and 
Ogata,— who also use a variational wave function simi- 
lar to ours, though with less variational freedom since it 
allows only c—f hybridization and antiferromagnetism. 
Indeed, we find that all additional variational parame- 
ters, e.g., /— / hopping and superconductivity, do not 
change appreciably the energy, hence the phase diagram. 
Moreover, the variational phase diagram bears similarity 
also to that one obtained by Martin, Bercx, and Assaad- 2 
by the dynamical cluster approximation, although in the 
latter case all transitions seem continuous. 



IV. CONCLUSIONS 

In conclusion, we have studied, by mean-held and vari- 
ational Monte Carlo techniques, the Kondo lattice model 
on a square lattice. The mean-held phase diagram is 
qualitatively but not quantitatively similar to the varia- 
tional Monte Carlo one. Restricting the analysis to the 
paramagnetic sector, we have found by variational Monte 
Carlo a large region of d-wave superconductivity, which 
however disappears from the phase digram once we al- 
low for antiferromagnetism. It is well possible that, if 
magnetic frustration is added besides that due to n c < 1 , 
superconductivity could emerge and intrude between the 
antiferromagnetic phase and the paramagnetic one, thus 
offering a possible explanation to what is observed in 
many heavy fermion compounds. 

The variational Monte Carlo phase diagram is prac- 
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tically the same as that obtained by Watanabe and 
Ogata— by a similar technique. In particular we also find 
that the onset of antiferromagnetism and the breakdown 
of Kondo effect are not simultaneous close to compensa- 
tion, n c = 1. Moreover, we find that the Kondo collapse 
always occurs via a first order phase transition and is ac- 
companied by a redistribution of low-temperature spec- 
tral weight at the chemical potential inside the Brillouin 



zone. 
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